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This  work  encompasses  two  geometric  evolution  equations  with  applications 
to  image  processing. 

First  we  show  the  existence  of  a  unique  solution  to  the  flow  of  an  H-system 
with  Dirichlet  boundary  condition  which  is  regular  up  until  the  first  time  of  energy 
concentration.  If  we  assume  the  solution  satisfies  a  certain  energy  inequality,  then 
the  solution  exists  for  all  time  and  is  smooth  except  at  finitely  many  singularities. 
The  behavior  of  the  solutions  at  these  singularities  is  also  discussed. 

Under  certain  conformality  conditions,  the  "H"  in  the  steady  state  version  of 
these  systems  is  in  fact  the  mean  curvature  of  the  solution.  This  equation  arose  from 
the  study  of  minimal  surfaces  with  prescribed  mean  curvature,  which  in  turn  came 
from  Plateau's  classic  "soap  bubble"  problem.  The  steady  state  system  was  studied 
in  detail  by  Hildebrandt,  Wendt,  Brezis  and  Coron,  et  al.  It  was  then  observed  that 
these  H-surfaces  (the  solutions  of  the  H-systems)  could  be  studied  in  more  generality 
if  they  were  viewed  as  asymptotic  solutions  to  the  flow  of  an  H-system.  Moreover, 

vi 
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due  to  their  structure  they  could  be  analyzed  using  similar  techniques  as  those  used 
to  study  harmonic  maps.  This  is  the  approach  that  we  take  in  our  analysis. 

Next  we  study  a  modified  version  of  the  mean  curvature  flow  equation  which 
can  be  used  as  a  model  for  image  restoration.  This  equation  is  interesting  since 
the  level  sets  of  the  solution  can  be  used  to  extract  details  from  an  image,  while 
the  mean  curvature  provides  for  anisotropic  diffusion  which  aids  in  noise  removal. 
We  have  introduced  a  control  factor  that  slows  the  evolution  of  the  level  sets  near 
significant  features  in  the  image.  We  discuss  precisely  how  the  model  creates  this 
behavior,  investigate  some  of  its  geometric  properties,  and  compute  a  bound  for  the 
extinction  time  of  the  level  sets  of  its  solution. 


CHAPTER  1 
INTRODUCTION 


1.1    Nonlinear  Evolution  Equations 

In  this  work,  we  study  two  different  nonlinear  evolution  equations.  The  equa- 
tions I  have  studied  here  are  not  only  nonlinear,  but  also  admit  additional  difficulties 
such  as  lack  of  boundedness,  lack  of  energy  inequality,  and  degeneracy.  While  these 
characteristics  make  them  more  appropriate  to  modern  applications,  they  also  present 
added  difficulties  into  the  analysis  of  the  equations. 

The  first  equation  I  studied  is  the  flow  of  an  H-system.  The  main  obstacles 
in  studying  this  system  are  the  unboundedness  of  it's  solutions  and  lack  of  an  energy 
inequality.  This  prevents  us  from  using  any  standard  techniques  to  analyze  the 
equation.  Instead  we  must  rely  on  very  delicate  estimates  which  require  extreme  care 
in  deriving. 

For  the  second  part  of  this  exposition  I  studied  a  partial  diflFerential  equation 
(PDE)  based  model  for  image  processing,  namely  the  mean  curvature  flow  equation. 
This  model  is  used  to  remove  noise  from  an  image  while  preserving  its  significant 
features.  In  addition,  the  level  sets  of  the  solutions  can  be  used  to  extract  features 
from  an  image.  The  model  is  extremely  efficient;  however,  it  is  also  nonlinear  and 
degenerate,  creating  further  complications  in  its  analysis  and  implementation. 

1.2    Problem  1:  The  Flow  of  H-Svstems 

One  of  the  classic  studies  in  partial  diff'erential  equations  is  that  of  harmonic 
maps.   In  1985,  Struwe  proved  the  existence  of  harmonic  maps  into  an  arbitrary 
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manifold  [20].  The  techniques  used  in  [20]  have  since  been  used  to  study  other  general 
second  order  evolution  equations  with  variational  structure.  Two  such  equations  are 
the  Landau-Lifshitz  equation  and  the  flow  of  an  H-System. 

In  the  following  work,  we  analyzed  the  flow  of  the  H-System.  This  flow  actually 
has  rich  geometric  interpretation.  If  the  solution  to  the  steady  state  H-system  satisfies 
certain  conformal  conditions,  then  " !!(«)"  represents  the  mean  curvature  of  u.  We 
prove  the  existence  of  a  unique  solution  which  is  regular  up  until  the  first  time  of 
energy  concentration  as  well  as  prove  global  existence  and  and  discuss  the  behavior 
of  the  singularities. 

1.3    Problem  2:  PDE-Based  Image  Processing 

Computer  vision  is  another  fascinating  subject  for  which  partial  diflferential 
equations  serve  as  extremely  eflFective  models.  These  models  work  in  a  systematic  way, 
performing  all  necessary  tasks  simultaneously.  They  remove  noise  while  retaining  and 
even  enhancing  significant  features.  The  model  we  study  here  arises  from  evolving 
level  sets  of  an  image  by  their  mean  curvature. 

This  model  is  based  on  the  work  of  Osher  and  Sethian  in  [18].  The  level 
sets  of  the  solution  to  this  equation  can  be  used  to  identify  important  features  in 
an  image.  Embedding  the  curves  as  level  sets  of  a  surface  allows  for  changes  in  the 
topology  of  the  curves  without  disrupting  the  evolution.  In  [9],  Evans  and  Spruck 
proved  various  geometric  properties  as  well  as  the  extinction  time  of  level  sets  to  the 
mean  curvature  flow  equation.  We  introduce  a  control  factor  into  this  equation  to 
inhibit  the  evolution  of  the  level  sets  near  the  boundaries  of  signiflcant  features  in 
an  image.  The  existence  of  a  viscosity  solution  to  this  modified  equation  has  been 
shown  in  [4].  We  study  how  this  control  factor  affects  the  geometric  properties  as 
well  as  the  extinction  time  of  the  level  sets  of  our  solutions. 


CHAPTER  2 
PRELIMINARIES 


2.1  Motivation 


The  equations  we  study  here  are  too  complex  to  arrive  at  explicit  solutions. 
Instead,  we  must  be  satisfied  in  studying  the  existence,  uniqueness,  regularity  and 
various  other  significant  properties  of  their  solutions.  The  natural  spaces  to  look 
for  and  study  these  solutions  in  are  Sobolov  and  Holder  spaces.  For  the  ease  of 
the  reader,  the  definitions  and  properties  of  these  spaces  that  will  be  needed  in  the 
following  chapter  are  included  here.  This  information  can  be  found  in  [11]  and  [17]. 

2.2    Sobolov  Spaces 
Throughout  this  chapter,  1^  c  M" .  For  convenience,  we  will  use  the  notation 

Jl',  =  nx[s,t]    and       =  fi*. 

As  usual,  C^{^)  will  denote  the  set  of  infinitely  differentiable  functions  with  compact 
support  on  fl. 

Definition  2.2.1  [Weak  Derivative]  Ifu,ve  Ll^{Q)  and  a  =  (ai,...,a„)  where 
the  ai's,  i  =  l,...,n,  are  non-negative  integers,  then  we  say  that  v  is  the  a^^-weak 
derivative  of  u,  written  v  =  D'^u  if 
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for  all  test  functions  ^  G  C^(f)). 

Fix  1  <  p  <  oo  and  let  k  and  /  be  non-negative  integers.  Define  the  Sobolov  space, 
W^(n)  as  follows: 


W!^{Sl)  =  {u  :        R  I  V|a|  <  k,  D^'u  exists  in  the  weak  sense  and  D^u  G  LP(f))} 


with  corresponding  norm 


E|a|<fceSS  SUp|D°u|,  p=QO. 


If  p  =  2,  we  write 


We  also  denote  by  H^{Q),  the  closure  oiC^{n)  in  H''{n).  In  fact,  one  can  show  that 


H^{n)  ^{ue  H''{n)\         =  Oondny  \a\<k-  1}. 


On  the  evolution  space  fi*  define  the  Sobolov  space 


W^'^inl)  ^{u:nl^R\  y\a\  <k,j<  I,   D^u,  Diu  e  L^(f]*)} 


with  corresponding  norm 


E\a\<k  (Jni  \Dy\'dxdtj  "  +  E     (4,  iDiul^dxdt)  \  l<p< 


oo 


E|a|<;fcess  sup|D^M|  +  J2j<i  ess  sup|D/u|, 


Hi 


p  =  oo 
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These  spaces  have  the  appropriate  structure  to  invoke  some  standard  results  from 
functional  analysis.  In  particular, 

Theorem  2.2.2  Wp{Q)  and  Wp'^{Ql)  with  the  above  defined  norms  are  reflexive 
Banach  spaces  for  any  1  <  p  <  oo. 

For  more  details,  see  [11]. 

2.3    Holder  Spaces 

Holder  Spaces  are  usually  the  appropriate  setting  for  studying  the  regularity 
of  the  solutions.  Although  it  is  not  always  possible  to  initially  find  the  existence  of  a 
solution  in  these  spaces,  one  can  find  a  solution  in  a  Sobolov  space  and  use  embedding 
theorems  to  conclude  that  the  solution  is  in  fact  Holder  continuous.  To  that  end, 
we  define  the  following.  Let  k  denote  any  positive  integer.  We  will  use  the  standard 
notation 

=  {m  :  1^  ^  M  I  D^u  exists  and  is  continuous  V  |7|  <  k}. 


Recall  that 


Mc'-in)  =  51  «ifxiD>|. 

|7|<fc 


For  0  <  a  <  1,  define  the  space 


with  corresponding  norm 


III  A  \u{x)  —  u(y)\ 
|w||co(n)  =  sup  '  ^  ^ 


x,yen  \x-y 


a 
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We  can  new  define  the  space 

C*+"(f})  4  {u  :  Q  ^  K  I  M  G  C\n)  and  L>>  G  C°(f7)V  I7I  =  k} 
with  corresponding  norm 


l«llc*+°(n)  -  ll^llc'cn)  +  W^^'^Wc-iU)- 


Finally,  on  the  evolution  space      we  can  define 

=  {u  :  Q*  ^  R  I  DlD'^u  e  C{nl)  \/2r  +  s<k  and 

D2u  e  c-{nl),D^^^u  e  Ct(fi^)V  I7I  =  A;} 

with  corresponding  norm 


\u\ 


2r+s<k      '  |7|=fc 


As  with  Sobolov  spaces,  Holder  spaces  also  afford  a  nice  structure. 

Theorem  2.3.1  C"=+"(f])  and  C*=+"'^(Q^)  with  the  above  defined  norms  are  Ba- 
nach  spaces. 

For  more  details,  see  [11]  and  [17]. 

2.4    Global  Approximations  by  Smooth  Functions 

If  one  begins  with  smooth  initial  and  boundary  data,  there  are  sufficient  tech- 
niques available  to  determining  the  existence  of  a  local  solution  to  the  nonlinear 
problems  we  study  here.  In  order  to  arrive  at  our  existence  results,  we  approximate 
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our  initial  and  boundary  conditions  with  smooth  functions  and  find  solutions  to  the 
system  using  this  smooth  data.  We  then  use  these  solutions  to  approximate  a  solution 
to  our  problem.  The  real  work  comes  in  finding  appropriate  uniform  bounds  on  our 
solutions  to  be  able  to  take  the  proper  limits.  Below  we  provide  the  approximation 
theorems  that  will  be  needed  once  our  estimates  are  suflftcient. 

Theorem  2.4.1  Assume  Q  is  bounded  and  u  E  Wp{Q)  for  some  1  <  p  <  oo.  Then 

there  exists      G  C°°(Q)  n  W!^{9)  such  that 

u""  in  W^{n). 

Theorem  2.4.2  Assume  Q  is  bounded,  dO,  is  C\  and  u  G  W^{^)  for  some  1  <  p  < 
oo.  Then  there  exists  u""  G  C°°{U)  D  W^{^1)  such  that 

u'^-^u    in  W^{n). 

See  [11]  for  proofs  of  the  above  theorems. 

2.5    Embedding  Theorems 

The  larger  the  space  in  which  we  search  for  solutions,  the  easier  it  will  be  to 
find  one.  Once  this  course  of  action  is  complete,  it  is  then  desirable  to  see  how  regular 
the  above  found  solution  can  be.  The  standard  way  to  do  this  is  using  embedding 
theorems.  In  the  following,  Bi  ^  B2  denotes  that  Bi  is  continuously  embedded  in 
B2.  As  in  the  previous  sections,     c  R". 

Theorem  2.5.1  1.  Suppose  1  <  p  <  00,  mp  <  n,  and  p  <  q  <  Then 


2.  Suppose  1  <  p  <  oo,  mp  =  n,  and  p  <  q  <  oo.  Then 
Furthermore,  if  p  —  1,  then 

3.  Suppose  \  <  p  <  oo,  mp  >  n,  and  0<A<m— ^  (<1)-  Then 

Remark  2.5.2  If  fl  is  bounded,  then  the  above  embeddings  are  compact. 
For  more  details,  see  [11]. 


CHAPTER  3 
THE  FLOW  OF  AN  H-SYSTEM 


3.1    Plateau's  Problem  and  H-Systems 

"Given  a  Jordan  curve,  F,  in      and  a  constant,  H,  find  a  surface,  u{xi,X2),  in 

of  mean  curvature  H  spanning  F." 

This  is  the  statement  of  Plateau's  classic  soap  bubble  problem  for  surfaces  with 
prescribed  mean  curvature.  One  might  recall  that  if  =  0,  then  u  represents  a 
minimal  surface. 

Analytically,  Plateau's  problem  can  be  stated  as  follows.  Let  M  C      be  a 
bounded  set  with  smooth  boundary.  Define 

M*  =  M  X  [s,  t]  and  M*  = 

A  map  u  e  C'^{M,R^)  satisfying 

Au  =  2H{u)u:,  Auy  (3.1) 

is  called  an  H-surface  (supported  by  M).  If  u  is  a  conformal  representation  of  a 
surface  5;  i.e.,  u  satisfies  the  relations 

then  H{u)  is  the  mean  curvature  of  S  at  u. 
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The  existence  of  /f-surfaces  under  various  boundary  conditions  have  been 
studied  in  detail  (see  [1],  [15],  [24],  [25]).  In  all  of  these  cases,  the  mean  curvature  was 
assumed  to  be  constant.  If  one  considers  the  flow  of  an  i/-System,  this  assumption 
can  be  generalized. 

3.2    Flow  of  an  H-Svstem 

In  the  rest  of  the  chapter,  we  will  study  the  existence  and  behavior  of  the 
singularities  for  the  heat  flow  of  the  /f-system  associated  with  (3.1), 

dtu  =  Au  —  2H{u)ux  A  Uy,  in 
'  w(0,a;)  =  uq{x),  in  M  (3-2) 

u{t,x)  =  x{x),  on{dMf 

where,  uo  G  H\M),  x  e  Hl{dM),  uo{x)\^3m)  =  x{x)  and  H  G  W^{R^). 

The  existence  of  a  global  regular  solution  to  (3.2)  has  been  shown  by  Rey  [19] 
assuming  uq  G  H^{M,R^)  n  L°°(M,M^)  and 

||-f^||L-(R3)||ito||L-{M)  <  1-  (3.3) 

The  assumption  of  small  initial  data  is  essential  for  the  proof,  since  this  leads  to 
an  energy  inequality  and  a  boundedness  property  for  the  solution  of  (3.2),  which  in 
general  does  not  hold. 
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Moreover,  Struwe  [21]  studied  the  free  boundary  problem  associated  with 

(3.2): 

dtu  =  Au  —  2H{u)ux  A  Uy,  in 
u{0,x)  =  uo{x),  in  M 

(3.4) 

u{t,x)  e  E,  a.e.  on  (dM)'^ 

dnu{t,x)  _L  T„(t,a;)E,  a.e.  on  (dM)'^ 

where  E  is  a  surface  diffeomorphic  to  the  sphere,  Tu(^t,x)'^  denotes  the  tangent  space 
to  E  at  u{t,x),  and  Ug  E  {u  e  H^{M)  :  u{dM)  C  S}.  Assuming  H  is  constant,  he 
showed  the  existence  of  a  local  regular  solution  to  (3.4).  The  free  boundary  condition 
complicates  matters;  however,  since  u\dM  ^  E  the  solutions  are  bounded,  a  property 
that  in  general  does  not  hold. 

Here  we  will  study  the  existence  and  the  behavior  of  the  singularities  of  (3.2). 
We  will  permit  i/  to  be  a  function  of  u  and  will  not  be  restricted  to  small  initial  data. 
We  will  show  the  existence  of  a  unique,  regular  solution  to  (3.2)  up  until  the  first 
time  of  an  energy  concentration.  In  general,  the  solution  does  not  satisfy  an  energy 
inequality;  however,  in  the  event  that  it  does,  we  can  show  the  existence  of  a  global 
solution  which  is  smooth  everywhere  except  possibly  at  finitely  many  singularities. 
We  will  conclude  by  discussing  the  behavior  of  these  singularities. 

The  main  difficulties  we  encounter  in  our  discussion  is  the  lack  of  the  energy 
inequality  and  the  unboundness  of  the  L°°-norm  of  the  solution.  These  issues  will  be 
addressed  throughout  the  remainder  of  the  paper. 


3.3  Notation 
Let  M  and       be  defined  as  in  section  2.2.  Define 

V{Ml)  =      e  C\[s,tlH\M))  :  \V^ul  \dM  e  L\MI)} 

where  the  derivatives  are  taken  in  the  distributional  sense. 
We  will  denote 

B^{x)  =  Br{x)  n  M  for  any  X  e  M 

where 

Br{x)  =  {x'  G  R^llx-x'l  <  R} 
If  z  =  {x,t)  e  M  X  [0,  oo),  then  for  t  >  i?^  define 

Qr{z)  =  Br{x)  x{t-  R\  t)  and  =  x  {t  -  R^,  t). 

We  will  denote  the  energy  of  u  as 

D{u;n)  =  1  [  \Vu\'^dx 

and  for  convenience,  D{u;  M)  =  D{u).  In  addition,  define  the  functional 

Eh{u)  =  D{u)  +  \  I  Q{u)ux  A  Uydx 

where 

H{S,  U2,  U3),J^     H{uu  S,  U3),  J      H{uu  U2,  s)]  . 
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Note  that  the  critical  points  of  Eh{u)  in  H^{M)  are  weak  solutions  of  the  steady 
state  of  (3.2). 


Finally,  the  quantity 


e{R)^    sup  D{u{t);B^{x))  (3.7) 


will  be  instrumental  throughout  the  remainder  of  the  paper. 

3.4    Local  Existence 
We  will  need  to  employ  the  following  Sobolov  type  inequality. 

Lemma  3.4.1  For  all  u  G  W^'^{M'^),  p  >  4,  we  have  that  Vm  €  C"'t(M^)  for 
a  =  {1  —  ^)  and  there  exists  Cp  >  0  such  that 


Theorem  3.4.2  Suppose  Uq  G  C^{M),  x  e  C°°(aM)  and  Uo\dM  =  X-  Then  there 
exists  an  5  >  0  and  a  unique  solution  u  G  C^+"'^+t(M'^)  of  (3.2)  on  M^. 

Proof:  For  any  p  >  4,  define  an  operator 

K  :  W^'\M'^)  ^  L^iM"^)  x  Wp'''{M)  x  Wp~'^{{dMf) 

by 

Ku  =  [dtu  -  Au  -  2H{u)u^  A  Uy,  u{x,  0),  u{x,  t)\^aM)T)- 
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The  differential  of  K  at  the  point  v  is 

DK{u)v  =  {dtV—Av-2H{u)[vxf\Uy+Ux/\Vy]—2DH{u)v-UxAUy,  v{x,  0),  v{x,t)\(^dM)'^)- 

By  [17] [IV,  Theorem  9.1],  there  exists  a  unique  solution  in  Wp'^{M'^)  to  the  corre- 
sponding parabolic  system 

dtv  —  Av  —  2H{u)[vx  A  Uy  +  Ux  /\  Vy]  —  2DH{u)v  ■  Ux  A  Uy  +  g  in 
^  v{x,0)  =  h{x)  in  M  (3-8) 

v{x,t)  =  k{x),  on  dM'^ 

for  any  {9,h,k)  e  IJ'{M'^)  x  Wp~'^{M)  x  Wp~^{{dMf).  Therefore,  DK(u)  is  an 
isomorphism. 

Moreover,,  evaluating  K  at  uq  yields 

Kuo  =  {go,uo,x), 

whereto  =  -Auo-2H{uo){uo)xA{uo)y  G  U'{M'^).  By  the  Inverse  Function  Theorem, 
there  exists  a  neighborhood  Ni  about  uo  in  W^'^{M'^)  and  a  neighborhood  iV2  about 
{go,uo,x)  in  D'{M'^)  x  Wp~'^M)  x  iyp~^((aMf )  such  that 


K  :Ni^N2 
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is  invertible.  Define  the  function 


95=  \ 


0  in 

f/o,    in  Mj. 


For  small  enough  5,  {gs,UQ,x)  ^  -^2,  so  there  exists  a  unique  u  e  Ni  e  Wp'^(M^) 
such  that  Ku  =  {gs,  uq,  x)'i  that  is,  u  satisfies 


dtu  —  Am  -  2H{u)ux  /\Uy  +  gs,  in 


'  u{0,x)  =  uo{x), 


u{t,x)  =  x{x), 


in  M 


on  {dMf 


(3.9) 


Hence,  there  exists  a  unique  solution,  u  G  Wp'^{M^)  to  (3.2).  By  lemma  3.4.1, 
Vu  G  C°'t(M'').    Finally,  applying  [17] [IV,  Theorem  5.2]  to  (3.9)  gives  us  that 


3.5    Some  A-priori  Estimates 


3.5.1    V(M^)  Estimates 

Lemma  3.5.1  For  any  smooth  bounded  domain  Q  C  E^,  and  function  (p  e  //^(Q) 


^  \<l>\'dx  <  cj^  \<l>?dx{j^ \V<j>\^dx  +  ^     101'^^^}  (3.10) 


with  a  constant  c  >  0  depending  only  on  the  shape  ofQ. 
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Lemma  3.5.1  can  be  found  in  [17]  [II,  Theorem  2.2  and  Remark  2.1].  As 
indicated  in  [20],  using  lemma  3.5.1  and  a  covering  argument  one  can  show  the 
following.  Refer  to  (3.5)  and  (3.6)  in  section  3.3  for  the  definitions  of  D{u)  and 
Eh{u). 

Lemma  3.5.2  There  exists  constants  c,Rq  >  0  such  that  for  any  T  <  oo,  any 
u  e  V{M'^),  any  R  G  (0,it!o] 

/    \Vu\^dxdt  <  cesssupD{u{t);  B^{x))- 

Jmt  {x,t)eM'r 


{[    \V\\^dxdt  +  R'^  [    \Vu\'^dxdt}.  (3.11) 

Jm'^  Jmt 

Moreover,  for  any  Xq  G  M,  any  R  e  {0,Ro],  any  u  G  V{M'^),  and  any 
function  rj  e  CQ°{Bji{xo))  depending  only  on  the  distance  \x-Xo\  and  non-increasing 
as  a  function  of  this  distance,  there  holds 

/    \Vu\^r]'^dxdt  <  cesssupD{u{t);  B^{xo))- 

Jm'^  0<t<T 

{       \V^u\Ydxdt  +  R-^  [    \Wu\Wdxdt}.  (3.12) 

Remark  3.5.3  We  will  henceforth  refer  to  Rq  determined  in  lemma  3.5.2. 
We  will  now  derive  L^-estimates  for  dtu  and  Vm. 

Lemma  3.5.4  Suppose  %  €  H\M),  x  G  Hl{dM),  uo{x)\^dM)  =  x{x)  and  H  G 
^LiR^)-  If  u  e  F(M^)  is  a  solution  to  (3.2),  then  D{u{t))  and  EH{u{t))  are 
absolutely  continuous  in  t  e  [0,T],  and 
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/    \dtu\^dxdt  +  E„{u{T))  =  Eh{uo).  (3.13) 
Proof:  Taking  the  derivative  of  D{u{t))  and  using  Greene's  identity,  we  obtain 


di 


D{u{t))=  [  Vu-Vdtudx=-[  Au-^tudxe  L^([0,r]). 


Using  this  fact  and  computing 


d_ 
di 


/  Q{u)ux  A  Uydx  =  3  /  H{u)ux  A  Uydtudx 


we  find  that 


^EH{u{t))=  [  [-Au  +  2H{u)uxAUy]-dtudx  =  -  [  \dtu\''dx  e  L'{[0,T]). 
Jm  Jm 

Integrating  this  last  equation  on  [0,  T]  gives  the  final  result.  □ 

We  will  need  to  employ  the  following  function.  Let  g  G  H'^{M)  satisfy 


Ai?  =  0,    in  M 

(3.14) 

5  =  X,      on  dM. 

By  the  theory  of  elliptic  equations,  there  exists  a  unique  solution  g  e  H^{M)  to 
(3.14)  such  that 
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(3.15) 


where  c  >  0  depends  only  on  M.  Recall  the  definition  of  e(i?)  given  in  (3.7)  in  section 
3.3. 

Lemma  3.5.5  Suppose  Uq  €  H\M),  x  e  Hl{dM),  Uo{x)\(^dM)  =  x{x)  and  H  G 
W^{R^).  Then  there  exist  constants  ci  and  ci  >  0,  depending  only  on  H  such  that 
for  any  solution  u  G  V(M^)  of  (3.2)  and  any  R  G  (0,i?o]  there  holds  the  estimate 


[    \V'u\'dxdt  <  c,D{uo)  +  cir(l  +  e,R~')  sup  D{u{t))  +  Cir||x|L. 

JmT  Q<t<T  H^{()M} 

provided  that  e{R)  <  ci. 

Proof:  Let  g  be  the  function  defined  in  (3.14).  Since  g  G  H'^{M)  and  u-g  e 
Hq{M),  then  by  the  Calderon-Zygmund  inequality  (see  [23]) 


/    \V^u\^dxdt<c{j     \u\^dxdt+  [    \Au\^dxdt  +  T\\g\\H2rM)). 
Using  (3.15)  and  Poincare's  Inequality,  we  can  estimate  the  above  by 

£^ \V^u\^dxdt  <  c( \Au\^dxdt  +  £^ \Vu\^dxdt  +  TM^^^^^^).  (3.16) 
To  find  a  bound  for  J^^r \Au\'^dxdt,  multiply  (3.2)  by  -Au  and  integrate  by  parts. 
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[   dt{^^^^)dxdt  +  [  \Au\^dxdt 
<  ||i/||L-  /  \Vu\^\/\u\dxdt 


By  lemmas  3.5.2  and  3.5.4,  the  above  inequality  becomes 


<-/     |Atxrd2;dt  +  c  /  \Vu\^dxdt. 


<  ce 


D{u{T))  -  D{uo)  +  \  f  \Au\^dxdt 

2  JM-r 


(R)  [    \V\\'^dxdt  +  ce{R)TR-'^  sup  D{u{t)).  (3.17) 

JM^  0<t<T 


Combining  (3.17)  and  (3.16)  we  have 

/    \V^u\^dx  <  cD{uo)  +  ce{R)  [  \V\\^dxdt 

+ce{R)TR-^  sup  D{u{t))  +  c  f    \Vu\^dxdt  +  cT\\x\\  3 

0<t<T  JmT  ^  ' 

<cD{uo)  +  cei        \V\\^dxdt +  cT {I +  eiR~^)  sup  D{u{t))  +  cT\\x\\  i 

JmT  0<t<T  H^{dM) 

Taking  €1  small  enough,  we  get 

f    iVM^dxdt  <  cD{uo)  +  cT{l  +  e,R-^)  sup  D{u{t))  +  cT\\x\\^^,,,,,.  □ 

J  M'^  0<t<T  «2(dM) 

Lemma  3.5.6  Suppose  uq  E  H\M),  x  G  Hl{dM),  uoix)\^aM)  =  x{^)  and  H  e 
^i(-R^)-  Then  there  exist  constants  C2  and     >  0,  depending  only  on  H,  such  that 
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for  all  solutions  u  G  V{M'^)  of  (3.2),  for  all  R  G  (0,i?o],  and  for  all  Xq  G  M  there 
holds  the  estimate: 

D{u{T);  B^{xo))  +  /  \W^u\^dxdt  <  2D(u(0);  B^j^ixo))  +  C2TR-'^c{R), 

provided  e{R)  <  €2- 

Proof:  Fix  Xq  G  M.  Let  77  G  Cq°{B2r{xo))  be  a  non-increasing  function  of 
the  distance  |a;  —  zqI  such  that  r/  =  1  on  Bii{xo)  and  [Vt^I  <  ;|  in  B2r{xo). 
Multiplying  (3.2)  by  -Au  ■  rj^  and  integrating  over       we  find  that 


/     dt{\\yu\^ri^)dxdt+  [    dtu  ■  Vu  ■  2r]Vridxdt  +  f  \Au\'^r]'^dxdt 


<  c 


<c  \Vu\'^\Au\r]'^dxdt 
/    \Vu\'^r]'^dxdt  +  1  [  \Au\Wdxdt. 


So  by  lemma  3.5.2, 


/    dt{\\Vufrf)dxdt+]-  f  \Au\Wdxdt 

<cf    \Vu\'ri^dxdt  +  S  [    \dtu\Ydxdt  +  c{5)  [  \Vu\''\Vr]fdxdt. 
Jmt  Jmt  Jmt 

<  ce{R)  [    \V\\Wdxdt  +  5  [    \dtu\Wdxdt  +  c{S)R~^  [  \Vu\''dxdt 

J'^'^  Jmt  J(Bil(x:^))T 


<ce{R)  [    \V\\Wdxdt  +  6  [    \dtu\Ydxdt  +  0(5)7 R-h(R).  (3.18) 
Jmt  Jmt 
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Since  u  is  a  solution  to  (3.2),  \dtu\  <  c{\V'^u\  +  iVup).  Using  this  fact  and  lemma 
3.5.2,  we  get  that 

/    \dtu\^rfdxdt<c  j     \V^u\^T]'^dxdt  +  c  \Vu\^rj'^dxdt 

JMi"  Jmt  Jmt 

<{c  +  e{R))  [    \V\\'^r]'^dxdt  +  cTR-h{R).  (3.19) 

Note  that  for  u  sufficiently  smooth,  integration  by  parts  yields 

/   [uxyl'^rfdxdt  =  -      UxU^yyrfdxdt  -  I  u^Uxy  •  2r]r]ydxdt 
Jm  Jm  Jm 

dxdt  +  /   UxUyy  ■  2r]r]xdxdt  —  /   UxUxy  •  2T]r}ydxdt 
Jm  Jm  Jm 

<2      \^u\Wdxdt  +  cl  \Vu\'^\\/T]\'^dxdt  +  l  [  \uxy\Wdxdt. 
J M  Jm  2  J 

Therefore, 

/  \uxy\Wdxdt<A  \Au\^T]^dxdt  +  c  [  \Vu\^\VT]\'^dxdt.  (3.20) 
Jm  Jm  Jm 

Hence,  approximating  it  e  V{M'^)  and  integrating  over  [0,r],  we  get  the 
estimate 

/     \V\\'^r]'^dxdt  <  c        \Au\'^r]'^dxdt  +  cTR-^e(R).  (3.21) 
Jmt  Jmt 

Using  (3.18),  (3.19),  (3.21)  and  lemma  3.5.5  we  find  that 
/    dt{\\Vu\W)dxdt+  [  \V'u\Wdxdt 

JmT  2  Jj^T 
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<  /    dt{-\Vu\Y)dxdt  +  I    \Au\Wdxdt  +  cTR-h{R) 

<  ce{R)  [    \V\\'^r)'^dxdt  +  c5  [    \dtu\'^T]'^dxdt  +  cTR-^e{R) 

<  c{e{R)  +  e{R)5)  [    \V\\Ydxdt  +  c{5)TR-h{R). 
Choosing  62  and  5  small  enough,  if  e(i?)  <  62  we  get 

f    dt{l\Vu\Y)dxdt+  [    \V\\^r]^dxdt  <  cTR-h{R). 

Therefore, 

D{u{T)-B^{xq))+  [    \V\\Wdxdt<  [  l\Vu{T)\'^rj^dx 

<  I  \\^u{<S)\Ydx  +  cTR-''e{R) 
Jm  2 

<  D(u(0);  B^j,{xo))  +  cTR-h{R).  □ 

3.5.2  Regularity 

Lemma  3.5.7  Let  u  G  C"+2,t+i(QM(^^))^  E  M  x  [0,oo),  be  a  solution  of  (3.2) 
with  uo  e  C°°(M),  X  e  C°°{dM)  and  uo\dM  =  X-  Then  there  exist  constants  63  >  0 
and  C  >  0  such  that  if  c{R)  <  63  then 

l|Vti|L-(QM(,^))  <  CR-^ 
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and  for  any  0  <  a  <  1, 


ll'^^llc'>-t(QM(,o))  <  C(l  +  \\x\\c->{dM)) 


(3.22) 


where  c  =  c{R,  a). 


Proof:  Let  ^  G  C, 


00/70)2 


be  a  cut-off  function  such  that  0  <  ^  <  1, 


^  =  1  on  Br{xo)  X  (^0  -  r^^o  +  ^^),  ^  =  0  outside  of  Br{xo)  x  {to  -  R'^,to  +  R/^), 
m\  <  for  |/|  <  2  and  |6|  < 

Set  [/  =  <.  Then 


Ut  =  AU-  f, 


'  U{0,x)  =  0, 


in  (zo) 


on  B^{xo)  x{t  =  to-  R^} 


(3.23) 


U{t,x)  =^{t,x)x{x),    ondM  X  (to-R^to) 


where  /  =  2Vm  •  -  uA{  -  2H{u)u^  A  u^^^.  Applying  [17][VII,Theorem  10.4]  and 
using  the  function  g  from  (3.14)  we  have  that 
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+  II|VmP^|Ilp(q^(zo))  +  llxllc-(aQ^(xo))) 
<c((i2-r)-ii??||Vu||^co(QM(,„)) 

+  (i?  -  r)-2  y|u  -  ^|Lp(qM(,„))  +  ||5|Lp(qM(,„)) 

<c  [{[{R  -  r)-'  +  {R-  r)-^]i2?||V«|L»(QM(,„))  +  rI\\Vu\\1^^q^^^^^^ 
+  [{R-r)-'Rt  +  l]\\x\\c-idQ-izo)))  ■  (3.24) 


By  lemma  3.4.1,  for  any  0  <  r  <  i?  and  any  p  >  4  and  a  =  (1  —  ^)  we  have 


<C||[/||^^2,1(qM(,^)) 

<C  ([(i?  -  r)-i  +  (i?  -  r)-^]i?F||V«||,oc(QM(,„)) 
+  i?p||Vn||2<»(^M(,„)) 

-  r)-2i2?  +  l]||xllc»(aQ-(.o)))  •  (3-25) 

where  C  =  C(a)  >  0.  Therefore,  it  only  remains  to  show  that 

ll^^llL-(Q-(.o))<Ci2-' 
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to  complete  the  lemma.  To  that  end,  choose  0  <  p  <  f  such  that 


-  pf  sup  |Vup  =  max  {(^  -  uf  sup  \Vu\^}  (3.26) 

2  Qp{zo)  0<<^<f      ^  Q<.(zo) 


Fix  zi  e  Qp(-2o)  such  that 


|Vwr(zi)  =  sup  |Vu|^  ^  e.  (3.27) 

Qp(zo) 


Then 


max  {(|  -  a)2  sup  |Vup}  =  (|  -  p)2e  (3.28) 


We  claim  that 


>  (3.29) 


By  way  of  contradiction,  suppose 


1  (f-^)' 


Define 


t;(a;,t)  =  ^(rri +e  ^^x,ti+e  4)  (3.31) 


and  denote 


Sr  =  gr(0)  n  {(a;,t)|(a;i  +  e'^^x,  t,  +  e'H)  e  Q^(zo)}. 


Then  v  e  C"+2't+i(5i)  and  it  satisfies 
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Vt  =  Av  —  2H{v)vx  A  Vy,    on  Si 

v{x,  t)  ^        +  e"^x),      if  xi  +  e~23;  G  dM  n  (^o)- 


Choosing  a-  =  p  +  e2         we  can  conclude  from  (3.28)  and  (3.30)  that 


sup  \Vv\^  <  sup       |Vup  <  e"^        sup        |Vwp  <  4.  (3.32) 

(i,t)e5i  (x,t)eQ^-i/2(2i)  (x,t)eo^^^_i/2(zo) 


In  addition,  (3.27)  implies  that 


\Vv\'{0)  =  e''\Vu\^{zi)  =  1. 


(3.33) 


In  order  to  arrive  at  our  final  contradiction,  we  only  have  left  to  show  that  there 
exists  a  constant  C  >  0  such  that 


'[  |Vt;| 


1<C  I  \Vv\^dxdt. 


(3.34) 


If  C  does  not  exist,  then  there  exists  a  sequence  {vi}  satisfying 


{vi)t  =  Avi  -  2H{vi){vi)x  A  {Vi)y,    in  Si 


(3.35) 


Vi{x,t)  =  x{xi  +  e  2x), 


ifxi  +  e-h  e  dMnQ^{zo), 


sup  \Vvi\^<4, 
{x,t)eSi 


(3.36) 
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\Vv,f{0)  =  1,  (3.37) 


and 


/   \Vvi\'^dxdt     0  as  i  ^  oo.  (3.38) 

Applying  (3.25)  to  Vi  (since  Vi  satisfies  (3.35)  and  (3.36))  and  taking  R  —  1  and 
r  =  |,  for  any  0  <  a  <  1  we  have 

W^^iWc-iis,)  ^  +  \\x\\c^{dM))- 

Therefore,  there  exists  a  subsequence  {fj^}  of  {vi}  and  a  function  v  such  that 

V^;i,  ^  Vt;  inC"''T(5i)  (3.39) 
where  0  <  a'  <  a.  Thus,  it  follows  from  (3.38)  that 

/    \Vv\^dxdt  =  0 
Jsi 


2 


which  implies  that  =  0  on  ^i.  However,  (3.37)  and  (3.39)  gives  us  that  |Vu|(0)  = 
1.  Since  v  6  ^''^(^i),  this  is  impossible  so  (3.34)  must  be  true.  Therefore,  if  (3.30) 
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holds,  then 


KC  Wvl'^dxdt 


[  |V^| 

JSi 

<  e-^C  [  \Vufdxdt 


<  e"^C     sup      /  \Vu\'^dxdt 

<  e~^C     sup      /  \Vu\'^dxdt 

te[to-R'^,to]  JB^ixi) 

Choosing  63  small  enough  leads  to  a  contradiction,  so  (3.29)  must  be  true.  Therefore, 
setting  (7  =  ^  in  (3.28)  and  using  (3.29)  gives  us  that 

f  f )    sup  |Vup  <       -  p")   <  4  (3.40) 


and  so 


l|Vu|L«,(QM(,„))  <  CR  1 


Finally  (3.25)  yields 


4 


where  C  =  C{R,a).  □ 
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3.5.3  Uniqueness 


Lemma  3.5.8  For  any  u  G  HliW')  it  holds  that  u  G  L'^{W')  and 


II'"IIZ,4(K2)   <  4||u||i2(K2)||VM||^2(K2). 


The  proof  of  lemma  3.5.8  can  be  found  in  [16]. 

Lemma  3.5.9  Suppose  Uq  G  H^{M),  x  £  H^{dM),  uq{x)\i^qm)  =  x{^)  (^nd  H  G 
W^{R^).  Further,  suppose  that  u,v  e  V{M'^)  are  weak  solutions  to  (3.2)  with  the 
same  initial  condition  u{0,x)  =  ^(0,  x)  and  the  same  boundary  condition  u{t,  x)\dM  = 
v{t,x)\dM-  Then  u  =  v  on  M. 

Proof:  Let  w  =  u  —  v.  Since  u  and  v  are  solutions  to  (3.2),  we  have 

\dtw  -  Aw\  <  \2H{u)ux  Auy  -  2H{v)vx  A  Vy\ 

<  2\H{u)  -  H{v)\\u,  A  Uy\  +  2\H{v)\\u,  A  {u  -  v)y\  +  2\H{v)\\{u  -  v)^  A  Vy\ 

<  \w\\Vuf  +  c\Vw\{\Vu\  +  \Vv\) 

Note  that  since  u,v  E  V{M'^),  then  for  a.e.  t  G  [0,r],  u,v  £  H'^{M)  ^ 
PF^'*(M).  If  we  multiply  the  above  inequality  by  \w\,  integrate  over  M,  and  apply 
lemma  3.5.8  we  get 


< 
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+c{f   \Vw\'')'2{[  (|Vu|^  +  |Vu|^)^(/  \w\^dxdt)'^ 

<c{      \w\'^dxdty^  +  c(  /    \Vw\'^dxdt)^{  /  \w\'^dxdt)* 
Jm^  Jm'  Jm* 

<  c||^i;||i2(M')||Vw||£,2(M*)  +  ^||Vt/;||i2(M,) 
<  c|k||i2(M')  +  l|Vu;||i2(M<) 

where  c  =  c{\\uo\\l2(m),  R,({R), \\u\\l^{m'), \\v\\l^{m'))-  After  integrating  over  [0,t]  for 
any  t  <T,  we  find  that 

/   ldt{\w\'^)dxdt  <c  [  [  \w\'^dxdt 

Since  w{0)  =  u(0)  -  v{0)  =  0,  the  above  inequality  becomes 

[  {\w{t)\^)dx  <  c  [  {[  \w\^dx)dt 
Jm  Jo  Jm 

By  Gronwall's  Inequality, 

/  {\wit)\^)dx<0-e^o^'''  =  0. 
Hence,  to  =  0  in  L'^{M^),  so  u;  =  0  a.e..  □ 

3.6    Existence  up  until  Time  of  Energv  Concentration 

Theorem  3.6.1  Suppose  uq  G  H\M),  uo\aM  =  X,  X  ^  Hl{dM)  and  H  e  W^{R^). 
Then  there  exists  a  unique  solution  u  E   n_V{M^)  of  (3.2),  defined  and  regular  on 
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M  X  (0,  T)  where  T  >  0  is  characterized  by  the  condition 

lim    sup  D{u{t)]B^{x))  >  e 

for  all  R  >  0,  with  a  constant  e  >  0  depending  only  on  H. 

Proof:  1.  Since  Uq  G  H^{M)  and  x  ^  Hi{dM),  there  exist  sequences  v]^  e 
C°°(M)  n  H\M)  and  x"^  e  C°°(aM)  n  Hl{dM)  satisfying  u'^{x)\sm  =  x'^ix)  such 
that 

wS"  ^  uo  strongly  in  F^(M)  (3.41) 

and 

X""     X  strongly  in      (dM).  (3.42) 

By  theorem  3.4.2,  there  exist  5m  >  0  and      G  F(M*™)  such  that  u""  solves  (3.2)  on 
M^""  with  initial  and  boundary  data  Mq*  and  x"^  respectively. 
Let 

e  =  min{ei,e2,e3} 

where  ci,  62,  and  63  are  determined  from  lemmas  3.5.5,  3.5.6  and  3.5.7  respectively. 
By  (3.41),  there  exists  some  R>  0  such  that  for  all  x  e  M, 

D{u^-B^^{x))<\. 
Then  by  lemma  3.5.6,  if  T  =  0{R^t)  then 
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sup  D{u"'{ty,B^{x))  <e.  (3.43) 

2.  Note  that  due  to  (3.43),  we  may  assume  that  5m  >  T.  The  reason  is 
as  follows.  U  Srn  <  T  for  some  m,  then  by  (3.43)  and  lemma  3.5.7,  we  have  that 
g  (j2+a,i+% ^j^5my  Therefore,  using  u^{x,  6m)  and  x"^{^)  as  initial  and  boundary 
data  for  (3.2),  we  can  apply  theorem  3.4.2  again  and  continuously  extend  the  solution 
u"^  to  a  larger  time  interval.  Since  we  may  keep  iterating  this  argument  as  long 
as  (3.43)  holds,  we  may  continue  the  solution  up  to  the  time  t  =  T.  Hence,  our 
assumption  is  safe. 

Using  lemmas  3.5.2,  3.5.4,  3.5.5,  and  3.5.6,  we  get  that 

[   ildtu^'l''  +  iV^u^^p  +  \Vu"'\')dxdt  +  sup  Diu"")  <  c{R).  (3.44) 

Jm'^  0<t<T 

Furthermore,  fixing  any  0  <  r  <  T,  for  any  Zq  6       and  Qr{xo)  C  M,  then  lemma 

4 

3.5.7  gives  us 

liv«»||^,,„.,,„„<c 

where  C  -  C{R,a).  For  any  compactly  embedded  M'  C  M,  a  standard  covering 
argument  results  in 

where  C  =  C{R,a,T,T). 

By  weak  compactness,  there  exists  a  subsequence  of  {it"*},  which  for  conve- 
nience we  will  still  denote  by  {u""},  and  a  function  u  e  V{M^)  such  that  as  m  ^  oo, 
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u""        weakly  in  y(M^).  (3.46) 
In  addition,  (3.45)  implies  that 


Vm"*  ^  Vu  uniformly  on  (M'^)  (3.47) 


for  any  M'  C  M. 

3.  We  claim  that  (3.46)  implies  that  u  is  a  weak  solution  of  (3.2)  on  V{M'^). 
In  fact,  by  (3.46),  as  m  — >^  oo  we  have  that 

V^"*  ^  V\  weakly  in  ^^(M^),  (3.48) 


dtu""     dtu  weakly  in  ^^(M^),  (3.49) 


and 


Vu"^  ^  Vu  weakly  in  L°°(0,  T;  L2(M))  (3.50) 
From  (3.48),  the  Sobolov  Embedding  Theorem  gives  us  that, 

VvT^Vu  strongly  in  ^^(0,  T;  L^(M))  (3.51) 
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and 


u'^^u  strongly  in  1^(0,  T;  L°°(M)).  (3.52) 

By  (3.48)  and  (3.49),  as  m  — )■  oo  we  have  that 

dtu""  -  Au""     dtu  -  Au  weakly  in  L'^{M^). 

So  it  only  remains  to  show  that  as  m  — >  oo, 

H{u'^)u'^  Au"^^  H{u)u^  A  Uy  weakly  in  L^{M'^).  (3.53) 

By  (3.52),  as  m  — >■  oo 

H{vJ^)      H{u)  strongly  in  1^(0, T;  L°°(M)).  (3.54) 

Furthermore,  (3.51)  and  (3.50)  imply  that 

<  ^     strongly  in  1^(0,  T;  L^{M)) 
and  weakly  inL°°(0,T;  L2(M))  (3.55) 

So  (3.53)  follows  directly  from  (3.54)  and  (3.55).  Therefore,  u  e  V{M'^)  is  a  solution 
of  (3.2)  in  the  sense  of  distribution  so  the  equation  holds  in  L^{M'^)  and  a.e.  in  fact 
it  is  a  solution  in  L'^{M'^)  and  hence  the  equation  holds  a.e.  Moreover,  from  (3.45) 
and  (3.47),  Vu  G  C°'t(M^^).  Then  by  [17][IV,Theorem  5.2]  u  €  C2+"'i+t(Mj)  is  a 
classical  solution  of  (3.2). 
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5.  This  solution  can  in  fact  be  extended  up  until  the  first  time  of  energy 
concentration.  This  is  due  to  the  fact  if  (3.6.1)  does  not  hold  for  any  x  E  M  a.t  time 
t  =  T,  then  (3.50)  holds.  Since  u  e  V{M^),  we  also  have  that  dtu  G  L'^{0,  T;  L'^{M)). 
Therefore,  u  €  C{[0,T];  H\M)).  In  particular,  u{T)  e  H\M).  Therefore,  the  above 
argument  guarantees  the  existence  of  a  solution  to  (3.2)  using  u{T)  and  the  new  initial 
data  and  we  can  continue  this  solution  to  a  larger  time  interval.  We  can  repeat  this 
argument  and  continue  the  solution  up  until  the  first  time  of  energy  concentration, 
that  is,  when  t  =  Ti.  Moreover,  the  solution  is  regular  on  M  x  (0,Ti]). 

Remark  3.6.2  If  x  ^  C^-^''{dM),  we  can  choose  Uq*  G  C°°(A/)  so  that  itQ^laM  —  X 
and  (3.41)  holds.  Then  by  lemma  3.5.7, 

ll^^llc°-t(Q|(zo))  ^       +  \\x\y+->{aM)) 

4 

where  C  =  C{R,a).  In  this  case,  we  can  conclude  that  u  e  C^'*""'^"'"f  (M^)  for  any 
0<T  <T  <T.  □ 

3.7    Global  Existence 

Theorem  3.7.1  Suppose  Uq  G  H^{M),  uo\dM  =  X,  X  ^  Hl{dM)  and  H  G  W^i(K3). 
In  addition,  suppose  that  for  any  0  <  ti  <  t2,  u  satisfies  the  energy  inequality, 

D{u{t2))  <  D{u{t,)).  (3.56) 

Then  there  exists  a  unique  solution  u  of  (3.2)  onMx  [0,oo),  which  is  smooth  on 
M  X  (0,  oo),  with  the  exception  of  finitely  many  singularities  {x[,Tk)  for  0  <  I  <  Lk 
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and  0  <  k  <  K ,  characterized  by  the  condition, 

limsupZ)(u(0;5fi  (4))  >  ^ 

for  all  R>  0,  with  e  depending  only  on  H. 

Proof:  Theorem  3.6.1  guarantees  the  existence  of  a  unique,  regular  solution  to 
(3.2)  on  M  X  (0,  Ti)  where  Ti  >  0  is  the  time  of  the  first  singularity.  By  the  additivity 
of  the  energy  and  (3.56),  there  can  only  be  finitely  many  singularities  at  that  time, 
{{x\,Ti)}\l^.  Therefore,  for  any  R  6  (0,i2o]  and  any  M'  C  M -iJ^^^B^{x\)  we  have 
that 


D(u(Ti),M')  <  limmfD(M(T),M') 


<  lmiinfD(u(T))  -  T,-^^,D{u{nB^{x\)) 


<  L>(u(0))  -  Lie. 


Letting  R^Q,  M'  ^  M,  one  can  conclude  that 


D{u{T{))  <  D{u{Q))  -  Lie.  (3.57) 

Since  D{u{T{))  <  D(u(0)),  we  have  that  u{T{)  e  H\M).  Therefore,  w(Ti) 
can  be  used  as  new  initial  data  for  (3.2)  and  another  application  of  theorem  3.6.1 
implies  that  the  solution,  u,  can  be  continued  to  a  larger  time  interval,  [0,72].  This 
process  can  be  repeated  at  each  time  u  has  a  singularity,  thus  obtaining  a  solution 
on  [0,oo).   Furthermore,  (3.56),  (3.57),  and  the  fact  that  D(u(0))  <  oo,  imply  that 
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there  are  only  finitely  many  times  {Tfcjj^i  at  which  u  can  attain  singularities.  Hence, 
u  can  only  have  finitely  many  singularities  on  [0,  oo). 

3.8    Behavior  of  Singularities 

Theorem  3.8.1  Let  u  :  M  x  (0,Ti]  be  a  solution  to  (3.2)  obtained  in  theorem 

3.6.1  with  smooth  boundary  data,  satisfying  (3.56).  Suppose  {xo,Ti),  xq  G  M,  is  a 
point  such  that 


for  all  R>  0  where  e  is  determined  in  theorem  3.4-2.  Then  the  following  holds. 

(i)  If  xq  G  M,  then  there  exist  sequences  tm      Ti,  x^  ^  Xq,  and  Rm  \  0  and 
a  non-constant  map  u  €  iZ/^g  fi  Cjg^iW')  such  that  as  m  ^  oo,  the  resettling  sequence 


\imsnpD{u{t)-B^{xo))  >  I 


(3.58) 


'^m(^)  —  ui^RjjiX  +  Xjji,  tjji^ 


converges  strongly  to  u{x)  in  Hj^^  n  C/^ 


7oc(^^)  (i'l^du  satisfies 


An  =  2H{u)ux  A  Uy  in  R^. 


Au  =  2H{u)ux  A  Uy  in 


< 


u{x)  =  x{xo) 


on  dRl 
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where  1^  =  {{x,x)\y  >  —a}. 

Proof:  First  we  will  prove  (ii).  Let  {xq,  Ti)  be  a  point  satisfying  (3.58)  where 
xq  e  dM  and  Ti  >  0.  By  the  finiteness  of  the  singular  set,  there  exists  a  5  >  0  such 
that 

u  e  C°°((Bf  (:ro)  x  [T,  -  5\T,])\{x,,T{)). 

Let  {Rm}  be  a  sequence  of  real  numbers  such  that  Rm  \  0.  By  lemma  3.5.7,  there 
exists  a  sequences  {xm}  and  {t^}  such  that  Xm     xq,  Tm     Ti  and 

sup        sup     /        \Du\'^{y,t)dy  =  e^  /  \Du\'^{x,Tra)dx.  (3.59) 

te[Ti-S^,Ti]  lesf  (xo)  JB^^ix)  JB^J^^m) 


Fix  0  <      <  2c2l{R)  where  C2  is  from  lemma  3.5.6.  Define 


Bm  =  {xe  W\Rn,X  +  Xme  (xq)} 


and 


Then 


Wm{x,  t)  =  u{RmX  +  X^,  Rl,t  +  T^). 


:  jB^  X  [-C^  0]  ^ 


satisfies  the  equation 


dtWm  =  -  2H{Wm){Wm)x  A  {Wm)y      OU  Bm  X  [-C^,  0] 

(3.60) 

Wm{x,t)=x{R^X  +  Xm)  ioTRmX  +  XmEdM 
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dxdt 


By  lemma  3.5.4,  we  have  that 

[  [ 

<  /   \dtu\'^dxdt  — >  0  as  m  ^  oo. 

So,  for  all  t  €  [-C2,0], 

/    \dtWm\^{x,t)dx — >■  0  as  m  ^  oo.  (3.61) 

JBm 

By  lemma3.5.6  and  (3.59),  for  any  f  6  [r^  -  C^R^,  nn], 


e  =  /  \Du\^{x,Tm)dx 


<  2  /  \Du\\x,  t)dx  +  C2{Tm  -  r)R-hiR) 

<  2  /  \Du\^{x,t)dx  +  -. 


Therefore, 


/  \Du\'^{x,t)dx>-. 

So  we  can  conclude  that, 


/    \DwJ'{x,t)dx=  [    \DuWx,Rl^t  +  Tm)dx 

JBm  JBm 
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Therefore,  for  all  t  e  [-C^,0], 


>[  \Du\'^{x,Rl,t  +  Tm)dx>^. 


/    \Dwm\'^{x,t)dx>^.  (3.62) 


Finally,  by  (3.59)  we  can  see  that 


sup    sup  /  \Dwm\'^{y,t)dy  =  (3.63) 

te[-c^,o]  xeBm  JBi{x)nBm 


sup  sup     /         \Du\'^{y,t)dy  ^ 

-C^Rl  ,Tm]  xeB^ (x„)  Jb^  (x) 


te[T„-c^Rl„Tm]  xeB^ixo)  Jb^^{x 

If  '^'^^(g^'^^)  ^  oo  then  fi^  ^  as  m  ^  oo.  Let  x  G  and  y  G  B^{x).  For 
m  is  sufficiently  large,  we  have  that  Rmy  +  G  M  so  by  lemma  3.5.7  and  (3.63)  we 
get  that 


sup  1 1  L>U;,ji  1 1 (Bi  {x)r>Bm )  <  c 
i6[-C2,0] 


which  gives  us 


sup    \\DwJ^\(.2+c,^2^<c.  (3.64) 

t6[-C2,0]  ^  ' 


Combining  (3.61),  (3.62)  and  (3.64)  we  can  find  r^^  e  (-C^O)  such  that  as 
m  — >  oo, 

/    \dtWm?{x,  T]jn)dx  — >  0  as  m  ->  OO  (3.65) 

J  Bm. 
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/    \Dwm\'^{x,r]m)dx>^  (3.66) 

and 

I \DWm{-,  Vm)  I  lcf+"(K2)  <  ^  (3.67) 

Therefore,  there  exists  a  function  u  :  — >  and  a  subsequence  of  Wjn{-,  ijm) 
such  that 

Wm{;  Vm)      u  strongly  in  Hl^  n  Cl^{^ ,  R^).  (3.68) 

Set  ^  =  r;^  on  both  sides  of  (3.60).  Letting  m  ^  oo  and  using  (3.65)  and  (3.68)  we 
can  conclude  that  u  satisfies  the  equation 

Au  =  2H{u)ux  A  Uy  in  R^.  (3.69) 

Moreover,  u  is  non-constant  by  (3.66)  and  (3.68)..  Therefore,  taking  tm  =  T^Vm  +  Tm 
we  have  that 

Vm{x)  =  Wm{x,  T]m)  =  u{RmX  +  Xm,  tm) 

is  the  sequence  desired  in  the  theorem. 

The  proof  of  xq  G  M  follows  an  analogous  argument  to  that  we  have  just 
presented,  so  it  will  be  omitted  for  brevity.  It  only  remains  to  show  (m)  case  2. 
Suppose  ^'^^(^--^^^      a  as  m  ^  oo.  Then 

B^^Rl  =  {{x^'),x^^))\x^')>-a}. 
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Note  that  on  the  set  {x^'^^  =  -a},  we  have  that  RmX  +  Xm^  Xq.  Also,  for  all  x  such 
that  RmX  +  Xm  e  dM,  Wm{x,t)  =  x{RmX  +  x^).  So  by  lemma  3.5.7  and  (3.63)  we 
can  conclude  that 


Therefore, 


sup     ||lOm||c2+°(Bi(x)nS„)  <  C. 

ie[-c2,o] 


sup     WWrnWcf+^iRl)  <  C- 


holds.  Employing  a  similar  argument  as  used  in  the  first  case,  we  can  obtain  a 
smooth,  non-constant  map  u  :      -^R^  and  a  subsequence  of  Wm{x,rim)  such  that 

Wm{-,Vm)^u  strongly  in  H^^^  n  C/o,(E^,  R^) 

and  u  satisfies  the  equation 


Am  =  2H{u)ux  A  Uy  in 
u{x)=x{xo)  on^R^.n 


CHAPTER  4 

MODIFIED  MEAN  CURVATURE  FLOW  EQUATION 


4.1    Evolution  of  Level  Sets  by  Mean  Curvature 

Consider  a  hypersurface,  c{t)  =  {xi{t),  ...,Xnit)),  evolving  in  the  direction 
opposite  its  unit  normal,  N,  with  speed  equal  to  its  mean  curvature,  k,  that  is, 

c'{t)  =  -kN. 

This  evolution  may  not  be  smooth  as  the  hypersurface  may  "pinch  off'  and  change 
topology  as  i  -)•  oo.  To  compensate  for  this,  we  can  embed  c{t)  into  a  surface,  u,  so 
that  c{t)  is  the  zero  level  set  of  u  at  time  t.  Let  w  :  R"  ->  [0,  oo)  be  defined  such  that 

c{t)  =  {xeW  I  u{x,t)  =  0}. 

Notice  that  on  this  level  set  of  u,  the  chain  rule  gives  us 

SILiWxi(a:i)f  +  Mt  =  0. 
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In  fact, 


Ut  =  -Vu  •  c'{t) 


--  -Vu  •  -kN 


=  Vu  •  div 


Wu\div 


Vu 


Vu 


Therefore,  the  evolution  of  the  level  sets  can  be  studied  by  considering  the  mean 
curvature  flow 


where  I{x)  is  a  smooth  surface  such  that  c(0)  =  {x  G  E"  |  I{x)  =  0}. 

This  equation  has  been  studied  in  detail  by  Evans  and  Spruck  in  [7]-[10]. 
They  proved  the  existence  of  a  unique,  smooth  solution  as  well  as  various  geometric 
properties  and  the  finite  extinction  time  of  the  level  sets  of  the  solution.  In  the 
following,  we  introduce  a  control  factor  into  the  right  hand  side  of  (4.1).  The  existence 
of  a  unique  viscosity  solution  to  our  modified  equation  was  proved  in  [4].  We  study 
some  geometric  properties  of  the  level  sets  of  this  solution  as  well  as  prove  that 
although  the  control  factor  may  slow  the  evolution  of  these  level  sets,  they  will  still 
eventually  shrink  to  a  point. 


< 


in  W  X  [0,  oo) 


(4.1) 


u{x,0)  =  I{x) 


in  R' 
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4.2    Applications  to  Image  Processing 

In  order  to  understand  the  applications  of  the  mean  curvature  flow  equation 
to  image  processing,  we  need  some  basic  definitions. 

Definition  4.2.1  Let  C  M"  and  TZ  C  R.  An  intensity  map  is  a  map  u  :  ^  ^  TZ 
where  u{x)  represents  the  intensity  of  the  image  at  the  point  x  ^  Q. 

For  2  or  3  dimensional  images,  we  usually  take  the  image  domain  to  be  f2  = 
[0, 1]^  or  [0, 1]^  respectively,  and  the  possible  range  of  intensity  values  to  be  7?.  =  [0, 1]. 
Most  PDE-based  image  processing  models  employ  Neumann  boundary  conditions. 
This  enables  us  to  consider  the  image  to  be  defined  on  all  of  R"  by  reflection. 

The  evolution  of  a  curve  as  described  in  section  4.1  can  be  used  to  identify  or 
extract  objects  from  an  image.  In  order  to  do  this,  the  curve  must  be  slowed  down  or 
halted  at  the  recognition  of  an  edge  (or  object  boundary)  in  the  image.  Furthermore, 
the  model  we  are  studying  here  can  also  be  used  to  retain  and  enhance  significant 
features  in  an  image  while  removing  noise  and  obstructions.  Therefore,  we  need  to 
formally  define  an  "object  boundary"  or  "edge".  Intuitively,  we  think  of  "edges"  as 
the  boundaries  of  the  significant  objects  in  an  image.  At  these  locations,  the  intensity 
of  the  image  should  change  significantly. 

Definition  4.2.2  If  u  is  the  associated  intensity  map  of  an  image,  then  an  edge  is 
any  place  where  |Vw|  achieves  a  local  maximum. 

Evans  and  Spruck  proved  in  [9]  that  the  level  sets  of  (4.1)  will  eventually 
shrink  to  a  point.  In  order  to  use  these  curves  to  extract  significant  features  we  need 
to  introduce  a  control  factor  that  inhibits  the  evolution  near  the  edges  of  the  image. 

To  that  end,  if  I{x)  is  an  image,  one  can  identify  its  edges  by  looking  for  large 
changes  in  intensity;  that  is,  local  maxima  of  |V/|.  We  do  not  necessarily  want  to 
assume  that  /  is  diff"erentiable;  in  fact,  one  would  hope  that  /  is  not  even  continuous 
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at  an  edge.  Therefore,  we  will  use  the  local  maxima  of  \VGc  *  /|  to  identify  edges, 
where  G„{x,t)  =  ;^^=fr  exp(^^)  is  the  Gaussian  filter.  This  enables  us  to  compute 
the  gradient  without  any  prior  knowledge  about  the  derivatives  of  /. 

Now  that  we  have  a  mechanism  for  locating  the  edges  of  an  image,  we  can 
construct  a  function  that  will  automatically  detect  these  edges.  Any  function,  g,  such 
that 

1,     away  from  an  edge 


0,     close  to  an  edge 


will  do  just  that.  Our  choice  is 


Introducing  g{x)  into  the  right  hand  side  of  the  mean  curvature  flow  equation 
(4.1)  provides  a  mechanism  for  slowing  or  stopping  the  evolution  of  the  level  sets 
as  they  approach  object  boundaries.  In  the  remainder  of  the  chapter,  we  study  the 
effect  of  the  control  factor,  g,  on  this  evolution. 


4.3    Anisotropic  Diffusion 

In  addition  to  the  evolution  of  the  level  sets  of  its  solutions,  the  mean  curvature 
flow  equation,  (4.1),  also  yields  a  special  type  of  difl"usion  which  is  very  useful  in  image 
restoration.  Let  ^  =  =  ^^^{-Uy,u^).  Notice  that  the  right  hand  side  of  (4.1) 
can  be  simplified  to 
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|v«MM^)  =  |v«M»(p^,^) 

"^■y^xx      '^^x^y^xy  ~^  ^x^yy 
"  |Vu|2 

Notice  that  the  second  order  term  is  simply  u^^,  the  second  directional  deriva- 
tive of  u  in  the  direction  of  ^  ||  Vu^.  Therefore,  (4.1)  performs  a  diffusion  only  in 
the  direction  perpendicular  to  Vu.  This  means  that  diffusion  will  only  be  performed 
in  the  direction  tangential  to  the  edges.  This  should  theoretically  prevent  blurring 
across  edges. 

Moreover,  the  control  factor  g  helps  regulate  the  amount  and  location  of  this 
diffusion.  Since 

\^u\div  {g{x)^^  =  g{x)\Vu\div  (j^)  +  ^9  ■  Vm, 

the  diffusion  term  is  inhibited  near  the  edges,  which  gives  us  another  mechanism  for 
preserving  significant  features.  In  addition,  Vg  •  Vu  is  a  hyperbolic  term  which  allows 
for  shocks.  These  shocks  encourage  discontinuities  at  the  edges,  which  is  of  course 
what  we  would  hope  for. 


I 
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4.4   Existence  and  Uniqueness  of  the  Weak  Solution 


Consider  the  problem 


ut  =  \Vu\div  [gix)^^      in    K"  x  [0,oo) 
u{x,  0)  =  I{x)     on  R" 


(4.3) 


where 


9{x)  = 


l  +  k\VG^*I\ 

In  [4],  Chen,  Vemuri  and  Wang  showed  that  there  exists  a  unique  bounded  viscosity 
solution  to  the  approximate  problem 


ul  =  {\Vu'\'^  +  e^)-2div(g(x)  ^^^^]      in  R"xfO,oo) 

u'{x,0)=I{x)       on  R". 


(4.4) 


provided  /  is  Lipschitz  continuous  on  R" .  Furthermore,  they  proved  that  as  e  ->  0, 


^  u  locally  uniformly  on  R"  x  [0,  oo) 


where  u  denotes  the  weak  solution  to  (4.3)  in  the  sense  of  viscosity  [5]. 

In  this  setting,  we  can  study  the  extinction  time  of  the  level  sets  of  u  as  well  as 
some  of  their  geometric  properties.  The  extinction  time  problem  can  be  reformulated 
as  finding  the  time  at  which  the  measure  of  the  level  sets  is  equal  to  zero.  Since  this 
is  always  the  case  if  we  use  Lebesgue  measure  to  determine  the  size  of  an  n  -  1- 
dimensional  object  in  R",  we  need  to  use  use  a  different  measure  for  defining  the 
extinction  of  the  level  sets. 
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4.5    Hausdorff  Measure 

We  are  accustomed  to  assigning  measures  to  n-dimensional  subsets  of  K" ,  the 
most  common  of  these  being  Lebesgue  measure.  However,  it  is  sometimes  convenient 
to  be  able  to  consider  measures  on  m-dimensional  subsets  of  E"  (m  <  n).  The 
following  information  can  be  found  in  Federer  [12]. 

In  1918,  F.  Hausdorff  introduced  an  m-dimensional  measure  on  K"  which  can 
be  defined  on  all  subsets  of  R" .  Moreover,  his  definition  is  consistent  with  the  usual 
notion  of  area  for  m-dimensional  submanifolds  of  R" .  The  definition  is  as  follows. 

For  any  subset  5  C  R" ,  define  the  diameter  of  S  to  be 

diam{S)  =  sup  \x  —  y\. 

x,yeS 

Let  am  be  the  Lebesgue  measure  of  the  closed  unit  ball  B(0, 1)  C  R"*.  For  small  S  >  0 
define  {Sj}j%^  to  be  a  countable  covering  of  A  such  that  Sj  C  R"  and  diam{Sj)  <  5 
for  j  =  1,2,3....  Then  for  any  A  C  ,  m  <  n,  define  its  m-dimensional  Hausdorff 
measure  to  be 

n''{A)^nm  inf  f  a^f^J^II^Y. 
For  any  a  e  ^  we  define  the  density  of  ^  at  a  to  be  the  quantity 

We  then  say  that  a  is  a  point  of  density  for  the  m-dimensional  Hausdorff  measure  if 
e"'(R'",a)  7^0. 

It  is  natural  to  restate  the  level  set  extinction  problem  as  finding  the  time  at 
which  there  are  no  points  of  density  on  the  zero  level  set  of  the  solution  to  (4.3). 
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4.6  Monotonicity 

In  order  to  further  discuss  the  geometric  properties  of  the  level  sets  of  the 
solution  to  (4.3)  we  first  need  the  following  monotonicity  estimate. 
Let     be  a  solution  of  (4.4). 

Theorem  4.6.1  (Version  1)  For  all  t  G  [0,T], 


Proof  :  Fix  S  >  0  and  define 


/        g{x)\Vu\dx  <  /  g{x)\VI\dx. 


0(x)  ^  e-''(i+l-l')^ 


for  all  X  eW.  Define 


Note  that 

\x\ 


^^'\t)  ^  [       g{x)ip^\Vuf  +  e^)Ux. 


(1  +  1X1^)2 


Using  the  notation 


i 


we  can  estimate 
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da; 


=   -  /        (f)'^div  g{x) 

-/ 

=  -/ 
-/ 

^  -/ 

VR"x 

/ 

^R"x 

/ 
/ 

JR 


+ 


+ 


R"x{f 
{t 
{t 
{t 
M"x{( 


(|Vn-|2  +  e2)i 


uldx 


2q{x)(b  r-u)dx 

2g{x)(j)V(f>-Vu'H'dx 
(P^H'YilVuf  +  e^y^dx 
g{x)(f)\H'f{\Vu'\^  +  e'')Ux 
g{x)\V(l)\''{\Vuf  +  e''f2dx 
g{x)cl)^{\Vuf  +  e^)'^dx 


Therefore, 


By  Gronwall's  inequality, 


[       g{x)(l>'{\W\'  +  £'y^dx<e'''  [  gix)cj)\\VI\'  +  e')Ux. 
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Letting  e,  (5  — >  0,  we  can  conclude  that 


/        g{x)\Vu\dx<  /  g{x)\VI\dx 


for  any  t  >  0.  D 

To  prove  the  that  the  level  sets  will  vanish  as  well  as  derive  a  bound  on  their 
extinction  time,  we  need  to  use  the  following  alternate  monotonicity  estimate. 
For  any  point  Zq  =  {xq,  to)  G  M"  x  [0,  oo),  define 

G,,{x,t)  ^  _i-exp(-—  -) 

V^47r(fo  -  t)  4(to  -  t) 

if  t  <  to  and      (x,  t)  =  0  otherwise. 
Theorem  4.6.2  (Version  2)  For  all  t  G  [0,T], 

/       g{x)\Vu\G,,dx<e'^'°'  [  g{x)\VI\G,,dx. 

Proof  :  For  convenience,  we  will  use  the  notation       as  defined  in  (4.5). 
Furthermore,  note  that 

and 

Define 

$^(0=  /  g{x){\Vuf+e'f^G,,dx. 


Then 


m'it)  =   /  g{x)  I  G.,dx 

+  /       g{x){\Vuf  +  e^)HG^o)tdx 

—    —  [        div  i  q(x)  r  I  ^^^G^n 

+  /       ^|io  -  Vnf  +  e^f^G.^dx 

-  /       six)       Jj"  u\VG,,dx 


Therefore, 


/ 


({t}  2(^0  -  t)' 

JK"x{i} 

+  /  '^\h-t\g{x){\Vu^f  +  e^Y^G,,dx 
+  /  ^^7^^MVu'WG,,dx 

jR-^x{t} 
JK"x{f}  ^ 


c{t} 

|2 


By  Gronwall's  inequality, 
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Letting  e  — >  0,  we  can  conclude  that 


/       9{x)\Vu\G,,dx<e('^'°y  [  g{x)\VI\G,,dxn 


4.7    Consequences  of  the  Monotonicity  Formula 


Consider  the  level  set  Ft  =  {x  E  R"  |  u{x,t)  —  0}.  In  this  section  we  derive 
estimates  on  the  size  of  Fj  for  any  t  >  0.  Moreover,  if  Fq  is  taken  to  be  a  non-smooth 
hypersurface  (for  example,  take  a  figure  eight),  then  Ff  may  develop  an  interior.  In 
this  last  case,  we  can  derive  bounds  for  the  size  of  dTt  for  any  t  >  0. 

Let  {^'fc}  be  a  sequence  of  smooth  functions  converging  locally  uniformly 
(away  from  s  =  0)  to 


Since  u  is  a  weak  solution  of  (4.3)  and  ^'jt  is  continuous,  v''  =  ^^(m)  is  a  weak  solution 
of  (4.3)  with  initial  condition  =  *jt(/)-  Therefore,  monotonicity  {Version  1)  holds 
for  v'';  that  is, 


1, 


s  >  0 


=  <     0,     5  =  0 


-1,    s  <  0 
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By  lower  semi-continuity  and  the  Lebesgue  Dominated  Convergence  Theorem, 


/        g{x)\Vv\dx   <   lim  inf  /  g{x)\Vv''\dx 
<    lim  /  g{x)\VvQ\dx 


lim  g{x)\VvQ\dx 
/  g{x)\Vvo\dx 


So  V  also  satisfies  a  monotonicity  property: 


/        g{x)\\/v\dx<  g(x)\\/vo\dx 


One  of  the  consequences  of  the  monotonicity  formula  is  a  certain  "  doubling  up" 

{l,s  >  0  0,s  >  0 

and  ^-(s)  =  <  .   Set  Vq  =  ^+(«o) 

0,s<0  I  l,s<0 

and  Vq  =  ^_(uo). 
Then 


/        \Vv\dx   <    -  [  g{x)\Vv\dx 
<  -  /  g{x)\Vvo\dx 

-/    g{x)\Vv+\dx+  j  g{x)\Vv^\d: 


<  -?^"-'(ronR")  + -?^"-i(ronE") 

C  C 
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where 


c  = 


1 


>  0. 


(4.6) 


1  +  ||-^||l°°(R") 


This  implies  that  Ft  has  finite  perimeter  for  any  t  >  0. 

In  addition,  if  Ft  is  the  closure  of  an  open  set  (this  may  happen  in  the  case 
that  Ff  develops  an  interior),  then  we  can  derive  a  bound  for  the  size  of  dTt,  namely 


Again,  consider  the  level  set  Fj  and  the  sequence  of  smooth  functions  {"^k} 
defined  in  the  previous  section.  As  in  the  last  section,  since  u  is  a  weak  solution 
of  (4.3)  and  is  continuous,  v'^  =  ^fc('u)  is  a  weak  solution  with  initial  condition 
Vq  =  ^fc(/).  Therefore,  monotonicity  {Version  2)  holds  for  u*;  that  is, 


We  will  now  show  that  even  with  the  control  factor  g{x),  although  the  motion 
of  the  level  sets  is  inhibited,  the  level  sets  will  still  vanish  in  a  finite  time  interval. 
This  is  primarily  due  to  the  fact  that  0  <  c  <  g{x)  <  1  where  c  is  defined  in  (4.6). 

First  we  show  that  if  only  a  small  portion  of  the  initial  level  set  lies  within  a 
ball,  then  there  exists  a  finite  time  at  which  the  evolved  level  set  no  longer  intersects 
a  smaller  concentric  ball.  This  result  is  analogous  to  Brakke's  "clearing  out"  lemma. 


4.8    Extinction  Times 
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Theorem  4,8.1  Suppose  Co  =  W'^iVo)  <  oo.  Given  p>  0,  there  exists  constants 
a,r]  >  0,  possibly  depending  on  Co,  such  that  if 

7{"-i(ron52p(0)<w"-^ 

then 

TtnBp{0  =  <Ii  fort^  [ap^2ap'']. 

Proof  :  For  simplicity  assume  p  =  1.  Let  t  €  [a,  2a]  where  a  >  0  is  still 
to  be  determined.  By  way  of  contradiction,  suppose  G  n  Bp{^)  is  a  point  of 
density  for  the  (n  —  1)  dimensional  Hausdorff  measure.  By  [13], 


lim  /        g{x)\Vv\G,^dx   >    lima  {VvlC^dx  (4.7) 

foVJjJnxld  to\t  yjRnx^tJ 

>   cialimr^-"  /       \Vv{-,t)\dx  >  C2  >  0. 


<{t}  ^0^^  ^R"x{t} 

r{xo) 

Then  by  monotonicity  and  (4.7)  we  have  that 


0   <   C2  <  lim  /  g{x)\Vv\G^odx 

to\t  yRnx{t} 

<    lime^*"*  /  g{x)\Vvo\G,,dx 
r  exp  (H^) 


<  dn^-^  +  n^-\ro)exp{--- 


v2 


V  47rQ:       \      e  8a 
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Choose  OL  =  a{Co,n)  >  0  and  rj  =  77(0;)  to  be  small  enough  so  that  the  right 
hand  side  is  less  than  C2.  Thus  the  result  holds. □ 

Finally  we  can  construct  a  bound  on  the  extinction  time  of  {rt}t>o- 

Theorem  4.8.2  Let  T  =  sup{t  >  0  |  7^  0}.  Then  there  exists  a  constant  C3  such 
thatO<  T  <  Ci?{"-nro)^. 

Proof  :  IfW-^iTo)  =  00  then  we  are  done.  Assume  that  0  <  n^'-^iVo)  <  00. 
Let  p  >  0  satisfy  ?^"~^(ro)  =  r]p^~^  where  r]  is  determined  in  theorem  4.8.1.  Then 
for  all  xo  G  M", 

n''-\TonB2,{xo))<r]p^-\ 
By  theorem  4.8.1,  this  implies  that 

Ft  n  Bp{xo)  =  0  for  ap^  <t<  2ap^. 

Since  this  holds  for  all     G  R",  we  have  that     =  0  if  i  >  ap^.  Therefore, 
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CHAPTER  5 
FURTHER  QUESTIONS 

5.1  H-systems  and  Related  Topics 

Many  interesting  problems  arise  from  the  analysis  of  the  flow  of  H-systems. 
Here  we  proved  the  existence  of  a  solution  which  is  regular  up  until  the  first  time 
of  energy  concentration.  Global  existence  was  shown  assuming  an  energy  inequality 
holds.  It  would  be  interesting  to  determine  if  there  are  weaker  sufficient  conditions 
that  provide  for  global  existence.  Another  issue  to  consider  is  the  nature  of  the  sin- 
gularities. For  harmonic  maps  it  has  been  shown  that  the  singularities  can  only  occur 
on  the  interior  of  the  region.  However,  it  is  still  not  known  whether  the  singularities 
of  the  Landau-Lifshitz  equation  or  those  of  the  flow  of  H-systems  occur  in  the  interior 
or  on  the  boundary.  Finally,  we  can  also  investigate  this  problem  in  higher  dimen- 
sions. There  has  been  extensive  study  of  solutions  of  harmonic  maps  from  E"  ->  ; 
however,  little  is  known  of  solutions  to  the  latter  equations  in  this  setting. 

5.2  PDE-Based  Image  Processing 

It  is  curious  that  using  (4.2)  as  the  control  factor  in  (4.3)  slows  the  evolution 
of  its  solutions,  yet  the  level  sets  still  shrink  to  a  point.  This  is  due  to  the  fact  that 

1  1 

~  1  +  A;|VG,*/|  -  l  +  ||/||L<»(Rn)  ^ 
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If  we  update  g  at  each  time  step;  that  is  define 

=  l  +  k\VG..u\  ^  °' 

then  (4.7)  may  fail  in  the  proof  of  theorem  4.8.1  and  it  is  not  clear  that  the  level  sets 
will  eventually  vanish.  It  would  be  interesting  to  see  if  in  fact  the  level  sets  do  not 
shrink  to  a  point  and  if  so,  determining  the  location  at  which  they  stabilize. 

There  is  also  much  analysis  to  be  done  on  other  existing  models.  Consider  the 

flow 

Ut  =  div{g{x)^~^^)  -  X{u  -  I) 

where 

=  i  +  ivk*/r 

Since  the  equation  is  in  divergence  form,  one  can  look  for  solutions  in  the  space  of 
functions  of  bounded  variation.  Furthermore,  one  can  also  see  that  this  equation  is 
equivalent  to  the  gradient  flow 

min  /  g{x)\Vu\dx. 
Jn 

In  this  latter  setting,  if  existence  of  a  solution  is  known  one  can  ask  about  the  long 
term  behavior  of  this  solution. 

Building  on  this  model,  consider  the  equation 

Ut  =  \^u\div{g{x)^^)  -  X\Vu\{u  -  I). 

This  equation  is  no  longer  in  divergence  form.  Due  to  this  fact,  only  the  existence  of  a 
viscosity  solution  has  been  shown  and  a  better  solution  is  not  possible.  Nonetheless, 
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we  can  still  investigate  geometric  properties  of  the  solution  as  we  have  done  in  this 
work. 
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